* tf-ii'-sS 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  of  This  PAGE  (Whan  D 


REPORT  DOCUMENTATION  PAGE 


RE  PCR  r NUMBER 

AIM  424 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2.  COVT  ACCESSION  NO.I  3-  RECIPIENT'S  CATALOG  NUMBER 


TITLE  (and  Subtitle} 

"The  Minimum  Energy  Movement  for  a Spring  Muscle 
Model  ^ _ 


TMORf»J  _ 

IJjJ John  MTHollerbach, 


5.  TYPE  OF  REPORT  & PERIOD  COVERED 


*.  PERFORMING  ORG.  REPORT  NUMBER 


(.  CONTRACT  OR  GRANT  NUMBERf»> 


9.  performing  organization  name  and  address 
Artificial  Intelligence  Laboratory 
5*»5  Technology  Square 

Cambridge,  Massachusetts  02139 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Advanced  Research  Projects  Agency 

f/i) 

1 400  Wi l son  B i vd 

Arlington,  Virginia  22209 

I 14.  MONITORING  AGENCY  NAME  8 APDRESS(7/  different  from  Controlling  Oltlca ) 

Office  of  Naval  Research 

J 

Information  Systems  ( ( Gr-J  ^ 

ftp.! 

Arlington,  Virginia  22217 

-7bJ 

(jS^  N000l4-77-C-jJ389"T (_ 


to.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  b WORK  UNIT  NUMBERS 


Se 


28 


S.  SECURITY  CLASS.  ( 

UNCLASSIFIED 


tS».  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (ol  (M>  Report; 


D i S^rtbut ion  of  this  d^iimpai-  is  unlimited. 

:J/£  v j 


17.  DISTRIBUTION  STATEMENT  (of  tha  abatrmct  antarad  In  Block  20,  It  dUfaraot  from  Raport) 


18.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Continu*  on  ravraa  old*  If  nacaaaa ry  ard  Idantlty  by  block  numb  or) 


D D 

?r^r‘  • 


energy  minimization 
muscle  model 
robotics 

motor  physiology 


modern  control  theory 


OP 


0 ABSTRACT  (Continue  on  reveree  tide  Hnecefew;  end  Idenllty  by  block  number) 

.There  are  many  ways  of  programming  an  actuator  or  effector  for  movement  be- 
tween the  same  two  points.  In  the  interest  of  efficiency  it  is  sometimes  de- 
sirable to  program  that  trajectory  which  requires  the  least  amount  of  energy. 
This  paper  considers  the  minimum  energy  movement  for  a spring-like  actuator 
abstracted  from  muscle  mechanics  and  energetics.  It  is  proved  that  for  this 
actuator  a bang-coast-bang  actuation  pattern  minimizes  the  energy  expenditure. 
For  some  parameter  values  this  pattern  is  modified  by  a singular  arc  at  the 
first  switching  point.  A surprising  limitation  on  the  duration  of 
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The  Minimum  Energy  Movement 
for  a Spring  Muscle  Model 


John  M.  Hollerbach 


ABSTRACT.  There  are  many  ways  of  programming  an  actuator  or  effector 
for  movement  between  th  same  two  points.  In  the  Interest  of  efficiency  it 
is  sometimes  desirable  to  program  that  trajectory  which  requires  the  least 
amount  of  energy.  This  paper  considers  the  minimum  energy  movement  for 
a spring-like  actuator  abstracted  from  muscle  mechanics  and  energetics.  It 
is  proved  that  for  this  actuator  a bang-coast-bang  actuation  pattern 
minimizes  the  energy  expenditure.  For  some  parameter  values  this  pattern 
is  modified  by  a singular  arc  at  the  first  switching  point.  A surprising 
limitation  on  the  duration  of  coast  is  demonstrated.  Some  relaxations  of  the 
restrictions  underlying  the  spring  model  are  shown  to  preserve  the  bang- 
coast-bang  solution. 
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Though  considerable  effort  has  been  expended  In  the  study  of  the 
human  motor  system,  the  execution  of  even  simple  movements  Is  not  well 
understood.  One  current  theory  holds  that  movements  are  memorized  In  terms 
of  final  position  tBizzl  et  al.].  The  organism  selects  length-tension 
curves  of  agonist  and  antagonist  by  adjusting  the  Innervation  levels  so 
that  the  intersection  of  those  curves  occurs  at  the  desired  position 
[Feldman] . 

This  process  is  illustrated  by  the  hypothetical  length-tension 
curves  of  agonist  and  antagonist  muscles  in  figure  1.  Suppose  the  system 
is  currently  at  length  Lo  under  innervation  rates  g(  for  the  agonist  and  nf 
for  the  antagonist.  If  the  Innervation  rate  of  the  agonlet  is  changed  to 
9a*  a different  agonist  length- tens  I on  curve  is  selected  and  the 
equilibrium  length  shifts  to  Lt.  Assuming  no  delay  in  tension  development 
and  ignoring  velocity  effects,  the  arrow  In  the  figure  Indicates  the 
tension  course.  There  is  an  Isometric  buildup  of  tension  from  P0  to  P2 
followed  by  an  .isotonic  decay  to  Pj,  where  the  tension  in  agonist  balances 
the  tension  in  antagonist.  The  final  position  theory  maintains  that  tho 
position  L|  can  be  reached  Independent  of  starting  position  merely  by 
sotting  rates  n|  and  g^.  This  theory  is  Interesting  from  a manipulation 
viewpoint  because  it  obviates  the  need  for  precise  trajectory  calculation. 

There  are  many  choices  of  agonist-antagonist  length-tension  curve 
pairs  that  have  L|  as  equilibrium  position.  One  choice  that  could  be 
expected  to  require  less  energy  Is  n2  and  g2,  which  minimize  the  ieometric 
tenelons.  More  generally.  It  Is  conceivable  that  some  complex  sequence  of 
Innervation  rates  (n|,g))  might  require  less  energy  than  a scheme  which 


?AQ£  3 


Figure  1.  The  equilibrium  point  of  the  intersecting 
length-tension  curves  of  agonist  (g  labels)  and 
antagonist  (n  labels)  shifts  from  LQ  to  1^  when  the 
firing  rate  of  the  agonist  is  raised  from  g^  to  g^ 
and  the  antagonist  rate  remains  at  n^. 


Figure  2.  Length-tension  curves  from  the  cat  soleus 
muscle  (Rack  and  Westbury) . 


7q  o a or  o o g 


••I  acts  the  final  length-tension  curvaa  1 mad  I Italy.  Tha  datarai  niton  of 
thli  optiaal  innarvation  pat  tarn  ii  tha  focui  of  thii  paper. 

I.  ICtueli  Muhanioa 

Tha  Invaitigatad  propart  las  of  auicla  praiant  a too  coaplicated 
view  for  analytic  treataent.  Tha  plan  In  ay  rmarch  hai  bean  to  ilapllfy 
tha  auicla  aachanlci  until  an  analytic  aolution  to  tha  optiaal  energy 
problaa  la  ponibla,  than  to  exsnlne  if  tha  natura  of  tha  oolution  la 
changed  by  adding  ioaa  of  tha  excluded  auicla  properties.  A full-blown 
foraulation  of  tha  problaa  would  require  nuaerical  Methods  for  solution, 
and  I intend  eventually  to  carry  out  this  analysis. 

The  length- tans ion  curvss  in  figurs  1 ars  theasslvss  a 
siapl  if  (cation  of  actual  length- tens  ion  curves  such  as  in  figure  2.  Tha 
siapllf  I cat  ion  arises  froa  extracting  only  the  linear  portions  at  short 
lengths.  Thera  is  soaa  question  as  to  what  portion  of  tha  length-tension 
curves  are  used  in  actual  aovsaent.  Soaa  authors  (Zierler,  Hill  1970,  Cook 
and  Stark]  aaintain  that  tha  active  portions  occur  near  Lo,  tha  length  at 
which  there  is  aaxlaua  Isoaetric  tension  and  which  Zierler  also  calls  the 
rest  or  natural  length  of  auscle.  Collins  at  al.  [19751  on  the  other  hand 
report  that  the  linear  portions  are  used  in  eye  aoveaent. 

For  the  present  we  accept  the  hypothetical  length-tension  curves  of 
figure  1.  The  curves  ars  also  assuaed  parallel,  as  reported  by  Collins  et 
a I . t Rack  and  Uestbury,  however,  report  a decrease  in  slope  with  firing 
rate  for  this  linear  portion  (figure  2).  The  curves  in  figure  1 lead  to  a 
nodal  of  auscle  as  a spring  with  variable  zero  setting.  The  slope  K of  the 
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curves  represent*  the  spring  constant,  and  the  variable  zero  setting  L, 
corresponds  to  the  selection  of  firing  rate.  The  force  exerted  by  a Muscle 
ie  thus  KCL-L,). 

An  iaportant  slept  If  Icetlon  of  ousels  properties  Is  to  exclude  the 
aeries  elastic  cooponsnt.  Ue  also  exclude  the  parallel  elastic  coaponent 
and  the  active  daaplng.  The  resulting  Muscle  Model  Is  presented  in  figure 
3.  The  aquation  of  Motion  for  the  spring  systao  of  figuro  3 1st 

a d^x  - - b dx  + k.Cx.-x)  - k„(x-xJ  Cl) 

dt*  « 

Define  a control  variable  X and  a state  variable  U as  below. 


. I x,  ) [ x ] 

X - l 1 • I ) 
C x,  J [ x ] 


[ u,  1 
U - I ) 
C u2  ) 


Setting  the  sees  a • 1,  the  state  variable  representation  of  the  spring 

eysten  1st 


. (0  11  [00] 

X - ( ] X ♦ I 1 U 

l *b  J ( kg  kR  ) 

flora  coapactly, 

X ■ A X + 0 U 
- f(X,U,t) 


(2) 


(3) 


II.  Mosols  Bnargatloa 

The  energy  E expended  during  aoveoent  equals  work  plue  heat.  The 
work  U May  be  subdivided  Into  conservative  work  perforaed  on  the  ness  a and 


Figure  3.  Simplified  muscle  model  with  equation  of  motion: 

m x = - b x + k (x  -x)  - k (x-x  ) 
g g n n 

where 

b ■ coefficient  of  passive  damping 
kg=  spring  constant  of  the  agonist 
^n=  spring  constant  of  the  antagonist 
m = mass 
x =■  position 

Xg=  variable  zero  setting  for  the  agonist 
x = variable  zero  setting  for  the  antagonist 
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nonconservat l ve  work  parforaad  on  the  viacoua  alaaant.  Tha  isoaetric  heat 
Qj  la  given  off  In  eatntaining  tha  auecle  at  a particular  tana  I on  P0.  Tha 
rate  of  energy  expenditure  la  thua 

a 

E • Po  v (power)  ♦ ct  P0  (aalntenance  heat  rate) 

■ (x4  + a)  (kf(u|-*|)  ♦ kn(X|-u2) ) (4) 


where  v la  velocity  end  a la  the  aalntenance  beat  coafflclant.  The  two 
force  terae  have  bean  aueeed  bacauae  each  contr  Ibutaa  to  energy  loea.  Ua 
have  excluded  the  ahortertlng  heat  bacauae  tha  active  daaping  was  also 
excluded,  and  bacauae  there  aay  be  a theoretical  relational) ip  between  the 
two  [Huxley,  Cap  I an) . The  tranaient  character)  at  lea  of  heat  production 
have  a loo  been  excluded. 


III.  The  Kular-Lagrange  Equations 

The  taak  now  ia  to  find  the  tlee  varying  control  U(t)  that 
• inieizea  the  energy  uaed  in  aovlng  between  two  points  in  a fixed  interval 
of  tiae.  Let  V represent  the  energy  conauaed  in  applying  the  control  U to 
yield  the  trajectory  X.  The  problea  of  ainiaizing  V ia  readily  approached 
by  technique  a of  aodern  control  theory.  Tha  fundaaental  equations  that  the 
optiaal  control  U(t)  aust  satisfy  are  derived  froa  a theorea  froa  the 
calculus  of  variations.  This  theorea  states  that  In  order  to  find  the  n- 
vector  XC t)  that  ainlalzas  V(X),  where 


V(X) 


*' 

L(X,  X,  t)  dt 


subject  to  the  constraint  relatione 

g,  (X,  X,  t)  - e 


(5) 

(6) 


I - 1 


a i n 
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then  X(t)  satiefiee  the  Euler  equations 

9V  (X.X.t)  - d_  3L*  (X.X.  t)  - 8 1-1 n (7) 

I dt  3k  | 


uhere 

L’(X,X,t>  - L(X, X,  t)  +Z  X|  (t)g,  (X, X, t)  (8) 

and  X|(t)'e  are  the  eultiplier  functions  [Schultz  and  falsa). 

Applying  this  the  ores  to  the  optisal  control  proto  lee,  the  state 

e 

equations  X - f(X,U, t)  represent  the  equality  constraints.  L represents 
the  rata  of  change  of  energy  E.  The  Haai  I ton  I an  H - L ♦ X f rapreeente 
(8),  uhere  X*  - tXj  Xfl.  By  applying  the  Euler  equation  first  for  X and 
then  for  If,  it  can  be  shown  that  the  alniaizing  U(t)  satisfies  the 
follouing  two  Euler-Laoranoe  equations  [Schultz  and  ftefsal. 


».-n, 

(9) 

"u  • e 

(18) 

IV.  The  Minimum  Principle 

Because  L is  linear  In  the  control  U,  there  uill  not  generally 
exist  a ninieue  energy  solution.  To  obtain  a realistic  solution, 
constraints  suet  be  placed  on  the  control.  The  solution  in  thie  case  uill 
Me  on  the  constraint  boundaries  (Bryson  and  Hoi.  Constraints  on  U, 
however,  make  it  impossible  to  differentiate  H with  respect  to  U. 

The  minimum  principle  of  Pontryagln  makes  it  possible  to  proceed 
from  thie  point.  Pontryagln  showed  that  even  if  the  control  is 
constrained,  one  etill  obtains  a minimal  solution  by  finding  the  u°  - 
u°(X,X.t)  to  minimize  the  Hamiltonian  H,  but  by  Inspection  rather  than  by 


i 
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differentiation.  After  finding  the  einieizing  u°,  one  fores  H°  • 
H(X,u°,A, t)  and  then  solves  the  following  two  equations  [Schultz  and 


Mai  sal. 


X - 3H° 

3JT 

A - - ai£ 

5 F 

There  are  tuo  natural  constraints  that  fall  on  the  control  U. 


First,  the  spring  cannot  push. 


U|  - *|  t 8 
x,  - uf  i 0 


Second,  springs  have  a aaxieue  tsnsion  that  they  can  exert.  Without  this 
constraint  the  solution  uould  involvs  an  Infinite  iapulse.  For  the  eoa**nt 
ue  assuae  the  eaxlnun  tension  is  constant  and  Independent  of  lengths 

(151 

*1  - S c2  (16) 

uhere  cj  and  c2  are  constants.  The  case  of  naxiaun  tension  varying  uith 
length  is  deferred  until  section  XI. 


V.  A Bang-Goast-Bang  Solatlon 

To  facilitate  inspection  of  the  Haei I ton i an,  ue  expand  H - XTf  ♦ L 
into  three  lines,  the  first  depending  on  ut,  the  second  on  uj,  and  the 
third  on  neither  control. 

H - kg(u|-x,)  (e+xj+X*)  (17) 

♦ («+Xf-Az> 

♦ Xj(X|-hX2) 
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To  atnlaize  H with  respect  to  U|,  we  observe  that  if  s+x^Xi  > ® then  H it 
ainlalzed  when  uj-k,.  If  *+xg+X2  < 0 than  H it  ainitlzed  ulth  U|«*|+c,. 
Stall  arty,  it  can  bo  shown  for  uj  that  when  a+x2-X2  < 0 the  ainiaizlng  u2 
Ilea  at  X|-c2i  otherwise  u2  is  at  x(.  Coabining  these  results,  one  finds  a 
bang-coast-bang  solution  to  the  ainiaue  energy  for  auscto  aoveaent. 

Case  1«  X2  < - («+xj) 

Then  U|-xt+C|,  Uj«xj. 

Case  2i  |X2|  i («+x2) 

Then  U|-X|,  Uj>X|. 

Case  3i  X2  > (#**2) 

Then  U|-X|,  u2-X|-c2. 


VI.  The  Solution  Equations 

Substituting  the  ainiaizlng  u°  into  H,  one  obtains  three  functions 
corresponding  to  the  three  cases. 

Case  1:  H°  - k,c,  («+x2+X2)  + x2(Xt-bX2) 

Com  Zt  H°  ■ *2(X|-bX2) 

Case  3i  0°  • k„c2(s+x2-X2)  ♦ x2(X|-bX2) 

The  differential  equation  (11)  and  its  solution  becoaes  for  the  three 

cases  i 


. (01)  ( 0 1 
X - I ) X ♦ t ] 
( 0 -b  ) t kc  1 
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»,<«)  - ♦ kcll  - 

b 

K,Ct)  - K,(t0)  ♦ l€xf(t0>  - Kt(tn  ♦ kc  (t-t») 
b T 

uh«r«  c-c,  and  k-k,  for  com  1i  c-8  for  case  2»  and  c— c2  and  k-k,  for 
caaa  3.  Tha  tlaa  t®  roproaanta  tho  starting  tlaa.  Tho  difforontial 
aquation  (12)  and  its  solution  arei 


. (001  (01 

A - ( 1 A - ( 1 

( -1  b 1 [ kc  1 


X|«t)  - A,(t0) 

A*(t)  - Af(t0)at,(t‘to)  ♦ Ai (tnlfkc  (1  - ,b,t’to), 

b 

uhara  c and  k hava  tha  saaa  aaaning  as  abova  sxcspt  c«+ct  for  casa  3. 
Sinca  A|(t)  is  constant,  it  appaars  bancs  forth  as  A|  ulthout  a tlaa 


(18) 

(19) 


(20) 

(21) 


VII.  Ths  Ixtrsnal  Varans  Singular  Solution 

It  Is  proven  In  append  Ik  A that  there  are  exactly  three  events  In 
the  extreaal  bang-coast-bang  solutions  an  acceleration  period,  a coast 
period,  and  a deceleration  period.  No  other  coabi nation  of  bangs  and 
coasts  is  ainiaizing.  However,  a nonextreaal  ainiaizing  solution  aag  arise 
froa  a singular  arc  at  the  switching  points.  The  Haailtonian  (17)  has  the 
curious  property  that  If  XzHa+xgl  then  the  corresponding  control  Bay  taka 
on  any  value  and  still  alniaize  H.  If  a control  can  be  found  to  aaintain 
Aflsr+Xfl  for  a finite  tiae  Interval,  then  a non-ex  trees  I solution  to 
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onergy  Minimization  Might  exist.  Thia  • I tuition  it  cal  lad  a singular  arc 
and  arlaaa  from  a performance  Index  linear  In  control  but  quadratic  In 
otata  tBrgoon  and  Ho).  To  maintain  Xt”l*+*tl  for  a finite  time  Interval, 
all  time  derivatives  of  the  two  switching  curves  Must  be  zerot 

d"(X»4«+x»)  • 0 nfc0  (22) 

dt" 

d"<X»-g-x»)  - 0 nZ0  (23) 

dt" 

Carrying  through  the  analysis  for  a singular  arc  at  ths  first  switching 
point  (22),  the  tine  varying  force  during  the  singular  arc  1st 


kghii-K,)  - X,  ♦ tw  ♦ 3bK1(t,)e2b(t"t,)  ♦ 3(ba+X()  (e2b(t"t,)-l)  (2*) 

2 


Unfortunately  no  sufficient  condition  has  yet  been  developed  to 
test  whether  a singular  arc  is  ninlnizing,  and  one  oust  coepare  values  of 
the  performance  Index  for  specific  parameter  values  for  the  singular  arc 
solution  versus  the  extremal  solution.  Depending  on  the  choice  of  Xi  the 
force  (24)  takes  one  of  the  three  forms  in  figure  4. 

Of  these  forms  only  4C  has  been  found  ninlnizing  for  sons  parameter 
combinations.  To  oearch  for  such  combinations,  a set  of  parameters  was 
Initially  deduced  from  Rack  and  Uestbury  (table  I).  The  elapsed  distance 
Xf  and  the  elapsed  tine  tf  are  variable  and  have  been  chosen  as  0.2  cm  and 
0.4  sec  respectively.  The  initial  and  final  velocities  are  assumed  zero. 
For  the  extremal  solution  (bang,  coast  .bang)  there  result  8 nonlinear 
equations  In  8 unknowns  from  (18) - (21 ) and  ths  initial  conditions  <Xf,tf). 
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For  the  singular  solution  (bang, 4C, coast, bang)  15  nonlinear  equations  In  15 
unknowns  result  froe  (18)  — (24)  and  the  initial  conditions.  The  equations 
were  solved  nueericatlg  by  Neuton-Raphson  and  gradient  Methods.  Individual 
paraeeters  were  varied  and  energies  of  eovseent  computed  froe  (5).  Solving 
(5).  the  energy  for  the  extreeal  solution  is 


E - kc«x,(t,)+3x,(t3))  + kca(at,+at3) 


(25) 


uhere  t(  is  the  switching  tiae  froa  acceleration  to  coast.  t3  is  the  tiae 
at  the  end  of  deceleration,  dtt  is  the  duration  of  acceleration,  dx((t|).  ia 
the  dletance  aoved  during  acceleration,  and  dx|(t3)  Is  the  distance  aoved 
during  deceleration.  For  the  singular  solution,  the  energy  is 


r*' 

4)  + / k,(U|-X|)  (X; 

A 


Z(t)-Hx)  dt  (26) 


E • kcOx|(t|)+ax|(t4>)  + kcaOt|+dt 

uhere  t(  is  the  switching  tiae  froa  acceleration  to  the  singular  arc  4C,  t3 
is  the  switching  tiae  froa  4C  to  coast,  and  t4  is  the  tiae  at  the  end  of 
deceleration.  The  force  kf(U|-X|)  is  given  by  (24),  while  the  velocity 
x2(t)  is 


f41  f4  . 2b(t-t|l  . . 2b(t-t,)  ,, 
x2(t)  • x2(t|)s  ♦ bfl^li Is  -1) 

The  energies  for  the  extreaal  versus  the  singular  solution  are 
coapared  in  tables  I la— gt  the  units  are  kg  ca/kg  ut.  In  table  lie 


(27) 


j 

i 
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Table  Ilflrfkiie.  b*4) 


st 

Singular 

Extremal 

0~06 

0.6112 

0.6161 

*0.1 

0.6647 

0.6681 

9.2 

0.7974 

0.7981 

9.25 

0.8633 

0.8631 

0.5 

1.1903 

1.1882 

1.0 

1.8393 

1.8382 

1.7 

2.748417 

2.748411 

1.8 

impossible 

2.8784 

th«  par a«e ter s k and  b are  varied  simultaneously  but  at  a fined  damping 
ratio  of  9.5.  In  tables  Hf-g  the  parameters  k and  b are  respectively  set 
at  16  and  4 rather  than  at  the  table  I values  where  the  extremal  solution 
la  minimizing  over  the  whole  range  of  x,  and  «.  The  Initial  values  from 
table  I are  starred  In  table  II.  A singular  solution  becomes  minimizing 
with  high  values  of  k,  b,  e,  and  tf,  and  with  low  values  of  at  and  x,.  As 
the  parameters  cause  the  coast  time  to  approach  zero  (higher  b and  xf> 
lower  k,  c,  and  tf),  the  singular  and  extremal  solutions  become  identical 
because  the  4C  portion  vanishes. 

For  the  extremal  solution  it  Is  proved  In  appendix  B that  there  is 
an  upper  limit  on  the  duration  of  coast.  It  Is  tempting  to  speculate  that 
for  longer  coast  durations  a singular  solution  becomes  minimizing,  but  the 
singular  solution  in  table  II  is  not  always  minimizing  under  these 
conditions.  Perhaps  a different  combination  of  bangs,  coasts,  and  singular 
arcs  would  then  be  minimizing,  but  this  remains  an  open  question.  Some 


combinations  can  be  proved  impossible,  such  as  (bang,  4C,  coast,  4A-C, 
bang) . 
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VIII.  Spring  Model  Relaxations:  X Dependencies 

A natural  question  Is  whether  the  Minima  energy  solution  is 
changed  by  incorporating  a More  realistic  auscle  Model.  For  those 
relaxatione  of  the  spring  Model  involving  only  X dependencies,  the  aneuer 
la  that  the  solution  reoalns  bang-cost t -bang.  The  reason  la  that  the 
Has  1 1 ton  I an  H reoalns  linear  in  the  control  U.  and  the  oinioization  of  M 
with  respect  to  U occurs  at  fixed  X.  Whether  the  solution  also  reoalns 
sees ! erst  i on-coast-decs  I erst  I on  nseds  to  be  deter  Mined  for  each  case. 

Relaxations  of  the  spring  Model  involving  X dependencies  include 
the  follouing. 

1.  Position  Halts  on  Tension 

For  real  Muscle  the  aaxiouo  isoaetric  tension  varies  with  position 
(figure  2).  This  Makes  c(  and  c2  into  functions  of  X|,  but  the  controls 
uill  still  fall  at  the  extreaes  wherever  they  are. 

2.  Velocity  Halts  on  Tension 

Actual  auscle  exhibits  a hyperbolic  force-velocity  relation.  If  P0 
is  the  isoaetric  tension,  then  the  Maxi  mum  force  P that  can  be  produced  for 
a velocity  v is  (Hill  19381  (see  figure  5)< 

P - P0  - viPn+s)  (28) 

v+b* 

The  terM  (P0+a) / Cv+b’ ) can  be  considered  the  coefficient  of  active  daaping. 
The  coefficient  a has  been  deteratned  as  .25  P0t  the  force  P then  becones 
P0-1.2SP0v/(v+b’).  The  literature  conflicts  on  the  value  of  active  daaping 
during  lengthening.  For  consistency  with  the  shortening  heat  (below) , it 


Tension 


>< 

•n 

•H 

O 

o 

rH 
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Figure  5. 


Tension  dependence  on  velocity  (Hill  1938) 


Figure  6 . Hypothetical  length-tension  curves  with 
the  property  that  at  any  given  length  the  slopes 
the  same  for  all  choices  of  u . 


are 
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ia  assuaed  tha  aaaa  aa  activa  daaping  during  shortaning. 

Aaaoclatad  ulth  tha  activa  daaping  la  an  extra  haat  expenditure 
abov*  tha  iaoaatric  haat  Am  to  ahortanlng.  Thla  abort  an  ing  haat  rata  ia 
[Hill  1964] t 

Q,  • .16  P0v  ♦ .18  Pv  (29) 

a 

Tha  Iaoaatric  haat  rata  Q|  roaaina  aP0,  but  tha  power  ia  now  Pv. 
Substituting  tha  aua  of  spring  for  css  for  tha  iaoaatric  tenaion  P0  and  (28) 
for  P,  tha  energy  rata  lot 

L - (k,(U|-x,)  ♦ k.luj-x,))  (a  - 8.89k,  ♦ 1.46b;)  (39) 

Xf+b* 

Slailarly  It  can  be  ohou  that  tha  aquation  of  notion  lot 

x,  ■ - b x,  + (1  - 1.2SxJ  (fc,(u.-x«)  - (^(xi-tif))  (31) 

XfJb* 

Uhan  thaao  taraa  are  coabinad  to  fora  tha  Haei  Itonian,  tha  control  la  aeon 
to  roaaln  linear.  Hanes  tha  solution  la  onca  again  bang-cooot-bang. 

»rfwg  Constant  Variation  with  Poaltlow 

Ota  way  of  bringing  tha  alapliflod  length-tension  curves  of  figure 
1 closer  to  those  of  figure  2 la  Illustrated  in  figure  6.  Tha  spring 
constant  k,  ear  las  with  position,  but  at  any  given  position  tha  constant  k, 
is  the  aaaa  for  all  controls  u,.  Under  these  conditions  tha  solution 
renal  ns  bang-coast -bang. 

4.  Parallel  and  Series  Elastic  Elaoanta 


i 


The  incorporation  of  those  a I soon  t a Into  tha  nodal  is  depicted  in 
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figure  7.  Since  the  parallel  elastic  eleeent  depend*  only  on  position,  it 
does  not  change  the  solution.  The  series  elastic  eleeents  and  the  active 
springs  My  toe  replaced  uith  equivalent  springs  uith  constants 
*** **V**s^  •h**  S’ "kjt+f (ka+k*) • This  eodification  also  has  no  effect 
on  the  solution. 

IX.  Spring  Modal  Relaxations!  U Dependencies 

In  figure  2 the  spring  constant  k„  is  seen  to  vary  uith  firing  rate 
at  any  fixed  poeitien.  The  linear  portions  of  these  length-tension  curves 
tdten  extended  seen  to  intersect  at  a common  point  (figure  8).  In  this 
circuestance  the  spring  constant  k„  is  VIuz-xq).  The  u2  teres  of  the 
Waal  I tonian  H becoeet 

H*  • Ta  (xi-u»)  («+x2-Xj) 
ur-xo 

- T0(«+x2-X2)  (xi-xn-1)  (32) 

U*-Xo 

If  s-fXf-Xx  > 0 then  H*  is  ainiaized  at  ug«X|t  otheruise  u2-x1-c2.  That  is 
to  say,  the  solution  for  U|  is  exactly  the  saae  as  in  section  V.  A siailar 
analyeie  holds  for  U|.  Thus  the  ainlaizing  pattern  is  also  bang-coast- 
bang. 

[ 


*TfcG6  A| 


Figure  7.  An  expanded  muscle  model  incorporating  series 

elastic  elements  k and  a parallel  elastic  element  k . 

s P 


Figure  8.  The  length-tension  curves  of  Rack  and  Westbury 

when  extended  meet  at  a point.  The  dependence  of  k on  u 

n 2 

can  be  characterized  by  T / (u  -x  ) . 
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II 

Appendix  A 

In  this  appendix  it  la  shown  there  are  only  tuo  bangs  and  one  coast 
in  the  extremal  eelutioni  one  acceleration,  followed  by  one  coast  period, 
terminated  by  one  deceleration.  No  other  combinations  of  coasts  and  bangs 
are  possible.  To  demonstrate  this  is  the  only  possible  combination,  it  is 
necessary  to  examine  the  switching  curves  and  their  time  derivatives. 

The  first  lemma  shows  that  once  the  control  has  passed  from 
acceleration  to  coast,  then  the  control  cannot  return  to  another 
acceleration  but  must  proceed  to  deceleration.  The  second  lemma  shows  that 
once  deceleration  has  started,  the  deceleration  must  continue  until  the  end 
of  the  movement.  This  proves  that  the  acceleration-coast-deceleration 
combination  Is  the  only  possible  one.  In  the  following  it  Is  presumed  that 
the  movement  starts  with  acceleration  in  the  positive  X|  direction.  Hence 
all  velocities  are  positive. 

Lemma  It  After  acceleration,  the  glide  period  cannot  double  back  to 
another  acceleration. 

Prooft  The  proof  of  this  lemma  proceeds  by  examining  the  time  derivative 
of  the  acceleration-coast  switching  curvs  (henceforth  referred  to  ae  the 
slope  of  the  switching  curve).  The  slope  of  this  curve  Is  initially 
positive  at  the  transition  from  acceleration  to  coast.  In  order  for 
another  acceleration  to  follow  the  coast  period,  this  slope  must  become 
negative,  leading  to  a contradiction. 

At  the  first  switching  time  tj  the  acceleration-coast  switching 


\L 


— 


curve  Is  zero 
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X*(t|)  ♦ « + *t<t|)  • 8 

After  the  acceleration  period,  the  coast  equations  erst 

% ...  * , btt-t|»  . M b(t-tj) 

X*t t)  ■ Xj(t,)s  ♦ ^(1  - a I 

b 


x2(t)  • - (X|(t|)-Hi)e 


— b 1 1— 1 1» 


* , b€t— t ,1  . M b(t-t,) 

J(|itl4e*K|it)  • Xjit|)s  ♦ ii.(l~e  I 

„ ..  . . — bit— 1|> 

♦ m - lX*lt,>+esls 


The  elope  of  this  switching  curve  isi 

ebU”t,,(bX2(t,)-Xl)  + blXl(t,)4«)e‘btt"t,)  (AS 

At  t«t|,  the  slope  of  the  switching  curve  1st 

2bXj(t,)-X|+b«  >0  (AS 

One  can  show  this  quantity  cannot  be  less  than  zero.  Next,  suppose  the 
coast  doubles  back  to  another  acceleration.  At  toes  point  the  slope  oust 
go  through  zero.  This  tine  t is  found  froa  (AS)  ast 

a2b(t”t,)  - b(X,(t,)+g>  > 1 (A) 

Xi“bX*(t|) 


Case  It  X|-bX2(t|)  > 0. 

Then  X2(t|)-Hi>0,  contradicting  (Al). 

Case  2i  Xi~bX2(t|)  < 0. 

Crossaultiplging  (A7)  and  collecting  teras, 
2bX«(t|)-X|-fbc<0.  contradicting  (AS). 


J 


. . « V. 
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Thus  after  accalaration,  tha  coast  period  aust  eventually  arrive  at  the 
deceleration  switching  point. 

Lemma  2\  The  Movement  is  locked  In  deceleration  until  the  end. 

Proofi  It  will  be  shown  that  if  deceleration  ever  switches  to  coast,  then 
the  elope  of  the  coast-deceleration  switching  curve  requires  an  imaedlata 
raturn  to  deceleration.  Hence  the  eoveaent  is  locked  In  deceleration  until 
the  end. 

Suppose  there  la  a time  t8  when  deceleration  switches  to  coast.  At 
this  point  the  coast-deceleration  switching  curve  is  zero, 

AjCta)  -a-x2(t3)  — 0 (A8) 


At  time  t*  the  elope  (A10)  la  bXt(t*)-&i+bxt(t8).  This  is  positive  since 
At<tS>>0'  *8(t8)>0,  and  A|<0  (lemma  3).  This  means  that  deceleration  would 
bounce  off  the  coast  boundary  and  immediately  continue  the  deceleration. 
Furthermore,  elnco  the  slope  is  positive,  the  deceleration  would  not 
Immediately  switch  back  to  coasting,  causing  chattering. 
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Lemma  3i  A|  < 8. 

Proof*  At  the  second  witching  point  tt  ue  have 

A*(t,)-«-)cfCtj)  - 0 (All) 

From  (A2)  and  (A3),  thle  becoeee 

A2(t,)eblt2"t,,+  - a + (^J(t,)4«)e*tolt*■t,,  - 0 (A12) 

b 

Rearranging, 

k'l-.6"*-''1.  . - »,(«,, (Al3, 
b 

Froa  (Al)  and  (A3)  ue  find  an  expression  for  Aa(t). 

Af(t|)  • - « - k€cj.(l-«”to  * * 1 ” ) (A14) 

b 

Substituting  into  (A13), 

Thus 

X,  - K.c,llVbltl~t»l)(.blt»-t'l«.-t<t»-<<l|  > b.ll«.bl>«-t'll  (Aiei 

1_.bct**tl) 

Since  the  nuaerator  is  positive  and  the  denominator  is  negative,  X(  ie 
negative. 

Taken  together,  these  leaaas  show  that  acceleration. passes  through  coast  to 
deceleration.  There  Is  no  posslbls  variation  in  this  scheae.  It  is  also 
possible  to  show  the  movement  cannot  start  bg  coasting  followed  by 
acceleration. 

n 

L — 


PAGE  26 


Appendix  B 


A surprising  lialtation  on  ths  vatus  of  tf-t|,  ths  duration  of  ths 
coasting  tins,  has  boon  found.  Ths  switching  cur  vs  during  accslsration  is< 


XfCt)4«4Kt(t)  - Xj(t0)sb(t"to,+ Xl±!Sa£i<l-«blt“t0,>  ♦«  <B1> 


Ths  slops  of  this  switching  curvs  1st 
b(t-t0) 


„ -htt-tg). 
♦ K,Ci  (1— s w ) 

b 


(bXjftoJ-Xi-kgCj)  ♦ AfCft 


-to(t-tg) 


(B2) 


At  ths  first  switching  tias  tj,  ths  switching  function  (Bl)  is  zsro. 
Rearranging  (Bl)  for  t»t|, 

b(t,-t0) 


(X2<t0)-(\i±^cL))  - - (\,+y:,)  - • - KTc,Cl-s~b(t|"t°)) 


(B3) 


b b 

Substituting  (B3)  into  (B2),  ths  slops  at  t(  1st 


-b(tj-tg).  . 

- 2k, c ] (l-«  ) - X|  - b« 


(B4) 


Substituting  for  X|  froa  (A16), 


..  -b(tf-t|)  b(t*-t|) 

K,c.  (1-a  ) (2+e  -s  ) + 2bg 

b(tj-tj) 


(BS) 


-1 


Sines  slops(t|)20  and  sines  ths  dsnoainator  is  positivs,  so  is  the 
nuasrator. 


■VlUV*,*rW>l2W*'‘r*,’-.,,,,r*,,l  . 2b.  t • IB6) 
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Bacauaa  2+a~b't2  t,*-ab,t2~t,,  ||  a dacraaalng  function  of  tx,  at  oom 
point  (B6)  bacoaao  zaro.  Solving  than  for 

ba  * V(l»Muc.ll-«~l’ltr<ol)>1  « ♦ 1 (B7> 

Aa  tj  Incraaaas,  tt  will  dacraaaa.  However,  tj  does  not  dacraaaa  anough  to 
offaat  tha  affact  of  tha  t2  incraaaa.  If  «-8,  (B7)  raducaa  to 

ab(t*"t,)  - 1 ♦ (B7a) 

Strangely,  In  thla  circuaatanca  t2-t ( dapanda  only  on  b. 
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